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We consider thermodynamic properties, e.g. specific heat, magnetic susceptibility, of alternating 
Heisenberg spin chains. Due to a hidden Ising symmetry these chains can be decomposed into a 
set of finite chain fragments. The problem of finding the thermodynamic quantities is effectively 
separated into two parts. First we deal with finite objects, secondly we can incorporate the frag- 
ments into a statistical ensemble. As functions of the coupling constants, the models exhibit special 
features in the thermodynamic quantities, e.g. the specific heat displays double peaks at low enough 
temperatures. These features stem from first order quantum phase transitions at zero temperature, 
which have been investigated in the first part of this work. 



I. INTRODUCTION 



This work is a continuation of our preceeding paper 
, which is devoted to the ground state properties of al- 
ternating Heisenberg spin chains. There we investigated 
Id systems described by Hamiltonians (model A) 



i/W = -J, J2 s(p) • {Tin) + <T(r2)) 

<p,r> 

-Jo '^iri)-Tir2), (1.1) 



<''l ,r2> 



or (model B) 



<p,r> 

<P1,P2> 

-Jo E cr(n) • T(r2). 



Two kinds of lattice sites, denoted by p and r, alternate 
within the chain. In both models, sites ri and r2 are 
occupied by nearest cr-spins (a = 1/2), which can be 
interpreted as forming a dumbbell configuration perpen- 
dicular to the chain direction, r denotes their common 
in-chain coordinate. In model A, p coordinates contain 
single s-spins, whereas in model B the p-sites are also oc- 
cupied by dumbbells of s-spins with coordinates pi and 
P2. A simple interpretation of model B is an alternating 
chain of orthogonal dumbbells. In this work we concen- 
trate on model A. However, the methods used below can 
be reformulated for model B as well. 

Two spins, cr(ri) and 'T(r2), are incorporated into the 
compound spin S(r) = <T(ri) -|- <T(r2), which is either 0, 
or 1. This reveals a hidden Ising symmetry of the orig- 
inal Heis enb erg model (|l . l|) . In fact, the Ji exchange 
term in (1.1) does not generate transitions between the 



total sp in s tates and 1 of any compound spin. Hamil- 
tonian (1.1) can be rewritten in a more suitable form as 



Hi+H^''\ where 



Hi = ~Ji E ^(p) ■ s(^) 

<p,r> 

and Ht^^-ljoJ2s'{r). 



(1.2) 
(1.3) 



Hq"^ counts the self-energy of a compound spin. 

We can use the following evident strategy: Any con- 
figuration of spins is characterized by intrinsic " defects" , 
i.e. r-sites, where the compound spin is zero. These "de- 
fects" , which are controlled by the Jo-terms, decompose 
the original chain into an ensemble of finite chain frag- 
ments, which are decoupled from each other. Their struc- 
ture can be defined as follows: A fragment of length k 
(fc > 1) is an alternating chain of fc + 1 spin-s sites and 
k spin-1 sites. Chain fragments are isolated from each 
other by a zero spin. 

It is convenient to measure all energies in units of 
Ji, the latter is supposed to be negative. Thus we set 
Ji - -1. 

In [|l| we have observed successive first order transi- 
tions governed by Jq at zero temperature. For model 
A with s = 1/2, it is a sequence (0) (1) (oo), 
where a periodicity element (fc) can be represented as 
(s,!)'^ , s,0. E.g., (0) is the periodical alternating struc- 
ture, where all r-sites are occupied by zero spins. For 
spins s = 3/2 and s = 2 the phase transition sequence 
becomes (0) ^ (1) ^ (2) ^ (3) (4) ^ (oo). The 
first two transitions, taking place at Jq^'^"^ and Jq^''^\ 

(2 3) 

are well isolated from each other, and from 

Jr '■ The 

latter appears to be very close to the values of Jq"^'^^ and 

„(4,oo) 
.7o 

A proper method to find out and classify all these 
transitions is based on linear programming theory. For 
our particular problem it prescribes to compare the re- 



1 



duced energies 



1 



1 



(cfe — kJo) of isolated chain frag- 



ments (s,l)'',s,0. €k is the ground state energy of the 
Hamihonian Hi (see ( |1.2| )) with open boundary condi- 
tions. It is convenient to introduce the following decom- 
position of Efe : 



Efc = fccoo + eo -I- ei„t(fc). 



(1.4) 



In (1.4) Coo is the energy per element (s, 1) of the per- 
fectly periodic spin structure, eo is the energy due to the 
open ends, and the remaining part, eint(fc), can be in- 
terpreted as the interaction between the chain fragment 
ends, which goes to zero at fc ^ oo. Thus, a succession of 
phase transitions is given by a broken line, which is con- 
cave upwards and envelops eint(^) from below Two 
typical functions eint(^) are shown in Figs, la and lb. 
In the former, eint(A) is a monotonic function, and thus 
the system passes through all intermediate phases from 
(0) to (oo), when jTo increases from large negative values. 
In figure lb, the enveloping function corresponds to a re- 
stricted number of transitions, which is typical for mixed 
Heisenberg chains , but with increasing values of s the 
minimum becomes very shallow. Note that if Hamilto- 
nian Hi takes the form of an Ising Hamiltonian, then 
eint(fc) = 0. In a spin-wave approximation for s > 3/2, 
this function belongs to the type shown in figure la. 



eint(k) 





FIG. 1. Two examples of possible shapes of eint vs 
k: (a) represents an infinite sequence of transitions, 
(0) ^ (1) ^ (2) ^ •■■ ^ (oo). In (b) a finite number of 
transitions, (0) (1) (2) (3) (oo), is realized. The 
dotted part of the line does not satisfy the condition "concave 
upwards" . 

Certainly, at non-zero temperature all phase transi- 
tions are smeared out. But thermodynamic quantities 
may exhibit a crucial dependence on temperature in 
vicinities of the critical j7o values. We start the next 
section with an example which will teach us some pecu- 
liarities of the thermodynamics of these systems. 

Besides analytic methods we have used numerical com- 
plete diagonalization of finite chain fragments. By uti- 
lizing global S"^ conservation we have been able to ob- 
tain all energy eigenvalues of the chain fragments with 
k — 1,...,6. The computations have been performed 



on an Ultra Enterprise 10000 computer manufactured by 
Sun Microsystems. This is a parallel computer with 40 
CPUs and 20 GBytes of shared memory. All programs 
have been implemented in CH — h. 



II. VICINITY OF THE (0) ^ (1) TRANSITION 

Let us consider the vicinity of the (0) (1) transition. 
For s = 1/2 and zero temperature the system under- 
goes the phase transition at j7o = — 2. At this value the 
ground state is manyfold degenerate: chain fragments of 
unit length, fc = 1 , are embedded into the (0) phase. This 
means that the distributions of "defects" (zero spins at 
r-sites) is subject to the following constraint: Two "non- 
defect" sites (spin-1) cannot be nearest neighbours on the 
r-sites. They must be separated by at least one "defect" . 
The partition function at zero temperature is the total 
number of all valid configurations. Let us assume that 
Zn counts all the configurations, which are realized be- 
tween the sites r = and r = n, but the sites r — and 
r = n are fixed at zero spin^. This yields a recurrence 
relation for the partition function Z: 



Zn — 2Z„_i + Zn-2 



(2.1) 



The two terms on the r.h.s. of ( p.l| ) correspond to the 
two possibilities for r = n — 1 to be either a " defect" or 
a "non-defect" site. If it is a "defect" site (first term), 
we count all configurations between the "defect" sites, 
r — and r = n—1. The factor 2 is due to the addi- 
tional spin-1/2 at site p = n — 1/2. This spin is enclosed 
by "defects", so it is free or paramagnetic. The second 
term on the r.h.s. of (^J) corresponds to the case, where 
r = n — 1 is a "non-defect" site. Since "non-defects" are 
not allowed to be nearest neighbours, r = n — 2 must 
contain a " defect" . Thus the configuration count of the 
remaining part is Z^-i. Evidently, the boundary condi- 
tions for (2.1) must be chosen as Zq = 1 and Zi — 2. 

Now we define an additional quantity which counts 
the number of " non-defects" between the spin-singlets on 
sites r = and r = n, summed over all allowed configu- 
rations. It satisfies the equation 



(2.2) 



The first two terms on the r.h .s. of (2.2) are similar to 
those on the r.h.s. of Eq. (2T). The third term counts 
how many times the " non-defect" at r = n — 1 appears in 
all valid configurations. For T>, the boundary conditions 
are = and Vi = 0. 



Let us discuss, how Eqs. (2J)-(2^) have to be modi- 
fied in case of non-zero temperatures. If all terms in Z„ 



^We label the sites of compound spins by integer numbers r. 
For p, half-integers are reserved. 
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are preceded by their statistical weights, we can use this 
quantity to calculate the partition function of a chain 
of length n with open boundary conditions. Instead of 
the set {efc} (cf. (1.4)), one should make use of the set 
of free energies {0/c}- Since we consider the vicinity of 
the (0) (1) transition (i.e. j7o — —2), only (j)o and 
(pi enter the calculation. All other configurations are 
suppressed at low temperatures (T <C 1). Introducing 
Boltzmann factors Wk = exp — ((/)fe — kJo)/T, and setting 

-2, we arrive at the following 



-Tin 2 and 



modification of Eqs. (|2.lD-(|2.2D 



2-2'^i 



wZ„ 



Wi, 



(2.3) 



This exhibits "half-a-gap" behavior, oc 1/y/w, when w ^ 
1, i.e., j7o > —2. This, and the analogous w-dependence 
of the specific heat can be interpreted in terms of an equi- 
librium chemical reaction, in which any spin singlet (1/2, 
1, 1/2) may transform into two paramagnetic spins 1/2 
and a "defect". 

In the vicinity of the (0) (1) transition, the results 
of this section are not only valid for s = 1/2, but also for 
other spins on the p-sites. Evident changes are: 

1. (pi — — (2s + l), and j7o varying around — (2s + l). 

2. (f)a = -Tln(2s + 1), resulting in 



I?„ = 2I?„_i + wiV„_2 + Z„_2), 
The solution of ( ^.3|) is 



(2.4) 



Z„ = c+Xl + c_ A!! , A± = 1 ± VTTw. 

For solving ( ^.4[ ), we try the ansatz 

V„ = {a+ + 72b+)Xl + (a_ + nb^)X'l. 

For the concentration of spin-1 sites in long chains we 
obtain 



lim 



nZ„ 



c+ 



(2.5) 



The relationship between and c+, which enter the 
leading terms of I ?„ a nd Z„ for n ^ 1, can be directly 
derived from Eq. (2.4): 



w 



w 



c+ XX+w 2{1 + w + yTTw)' 

The physical quantities, such as the specific heat 
and the entropy, can straightforwardly be calculated as 
derivatives of the free energy. In the thermodynamic 
limit, i.e. n oo, the free energy per compound spin 
is simply 



-Tin A I = - 



2 + Jo 



T 



w 



if (1) is favorable {w 3> 1). Otherwise, for ti; ^ 1 we 
arrive at the expression 



-rin2 



Tw 



The magnetic susceptibility reflects the groundstate in- 
vestigations of our former paper. In fact, as far as chain 
fragments of length k = 1 are in the spin-singlet state, 
only "isolated" spin-1/2 sites, associated with fc = 0, con- 
tribute a Curie-like susceptibility. It can be expressed as 

a contribution of individual "isolated" spins, multi- 
plied by their concentration 



and 



A± = (s + 1/2) ± v/(s + l/2)2 + w, 



w/2 



1-26+/C+ = (1 + Vl + w)/{l + W + Vl + w). 



c+ (s + 1/2)2 + W+ (5 + 1/2)7(1+172)2+^ 

3. A group of spins (s, 1, s), whose total spin value at 
low temperatures is practically (2s — 1), is param- 
agnetic too, as well as an "isolated" spin s. 

The Curie-like susceptibility is straightforwardly calcu- 
lated as 

^(2s(2,s-l)6+/c+ + .s(s + l)(l - 26+/C+) ). 



III. GENERAL CONSIDERATION 

For s = l/2, quantum fluctuations are efficient enough 
to "isolate" the (1) ^ (oo) transition from (0) ^ (1). 
In Ref. [0 these zero-temperature transitions were esti- 
mated as J'q^'°°-' — —0.910 and Jq^'^'' — —2, respectively. 
Thus, at low temperatures we can investigate the regions 
around J'q^'°°^ and Jq^'^^ separately. In Section II the 
subject of interest is the competition of "defect" and 
"non-defect" sites, provided two "non-defects" cannot 
be nearest neighbours, if Jo is around Jq*''^^ On the 
other hand, two "defects" cannot be nearest neighbours 
in the second critical range around Jg^'°°'' at low temper- 
atures. Instead of dealing with "defect" and "non-defect" 
objects, let us use the convention of Ref. ||l|. For conve- 
nience, chain fragments have been defined as follows: A 
chain fragment of length k formally includes a spin-0 site 
from its right, so it can be represented as (l/2,l)'^(l/2,0). 
Conventionally, the nearest spin from the left of any chain 
fragment is also 0, but this spin is incorporated into the 
fragment which lies on the l.h.s. of this site. In this clas- 
sification, a chain fragment of zero length is (1/2,0). 

Now we can reformulate one of the statements men- 
tioned above. In the second "critical" range around 
Jq^'°°^ zero-length chain fragments are incompetitive and 
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may be neglected at low temperatures. However, we shall 
see from the specific heat calculation that at intermedi- 
ate temperatures, e.g. T = 0.3, and intermediate JTq 
between i/q"'^^ and J^^^'°°\ not only chain fragments 
of unit length dictate thermodynamic properties, but 
longer chain fragments and zero-length fragments also 
contribute essentially. Therefore we have to take into 
account chain fragments of any length, and no special 
restrictions on the values of j7o and temperature will be 
imposed. 

The free energy of a chain fragment of length k can be 
written as a generalization of Eq. (1.4): 



0fc = fc/oo + /o + /int(fc), 



(3.1) 



where foo, fa and fmt are temperature-dependent. In this 
section, the configurational part of the free energy will be 
determined by making use of recursive relations similar 
to (P)-(|2l). 



It is convenient to take the global configuration (oo) as 
the "vacuum" state. A zero spin within this background 
i s ca lled a "hole" . If we ignore the interaction term in 
(34), i.e. /int(fc), then any "hole" costs the free energy 
(j7o — foc,{T)) + fo{T). In fact, the free energy of a very 
long chain, nf^{T) + fo{T), becomes (n - l)/oo(T) + 
2/o(T) -I- Jo if a "hole" is inserted. For two "holes" we 
obtain {n - 2)foo{T) + 2>fo{T) + 2 Jo, and so on. How- 
ever, this consideration is no more valid if two "holes" 
occupy nearest neighbour sites on the r-sublattice. Such 
two "holes" give rise to a zero length chain fragment. 
In this case we obtain (n - 2)f^{T) + 2/o(T) -I- 2Jo, 
which can be subdivided into {n - l)/oo (T) + 2/o(T) -t- Jo 
and — foo{T)Jo. The latter should be interpreted as the 
free energy of the zero length chain fragment. Only the 
"holes" which are separated by a chain fragment of non- 
zero length and which do not have other "holes" between 
them, interact via /int(fc;T), where fc > 1 is the chain 
fragment length, or the number of spin-1 sites between 
these nearest holes. We denote the statistical weight of a 
chain fragment of length k hy Wk- Then starting with 

Wo = exp [-(Jo - foo)/T] and 



= exp [-(Jo - /oo + /o + /int(l))/r] , 

we obtain 

zwfc =«;iexp[(/int(l)-/int(fc))/T], k>\ (3.2) 

for longer chain fragments. 

Recurrence relations f or t he partition function are ev- 
ident generalizations of (2.3): 



Z„+i = 2wo2„ + WiZn-l + W2Zn-2 

H h WkZn-k H \- WnZa, n>l. (3.3) 

Each Zk counts all possible spin configurations with cor- 
responding statistical weights between r = and r = k, 



while fixing the boundary compound spins at r = and 
r = fc at zero. The prefactor 2 of w^Zn on the r.h.s. of 
Eq. (3.3) is due to the spin-1/2 at site p — n + 1/2. The 
boundary conditions for the set of partition functions are 
Zq = 1 and Zi = 2wq. The latter reflects the existence 
of a free spin-1/2 between two spin-0 sites. 

The lower index in 2?^™^ has the same meaning as in 
Zk , whereas the upper index is related to the chain frag- 
ment length. I'^™'* measures how often chain fragments 
of length m occur between the sites r = and r — k. Of 
course, any spin configuration in 2?^™'' picks up a corre- 
sponding statistical weight. The recurrence relations for 
shortest chain fragments, fc = and 1, have a structure 
which is similar to that of Eq. (^.3[): 



V. 



(0) 



n+l 



2wQZn + 2woVl°^ 



(0) 



n > 1. 



n > 1. 



(3.4) 



(3.5) 



The first terms on the r.h.s. of ( |3.4|) and (3^) are the 
contributions of the spin-0 and spin-1 sites at r — n, 
respectively. The boundary condition, which should be 
imposed on reads V^^^ — 0. For 'D^^\ we can set 



V^Q^ = and V\ 



(1) 



0. 



Generalization of the recurrence relations and bound- 
ary conditions to arbitrary chain fragment lengths m is 
also evident: 



• • • + "'fe^«-fc + 



• • • + Wni^Q , n > m, 
Pi"' =0 for n<m. 



(3.6) 



We define the concentration Xm of chain fragments of 
length m as the ratio of the expectation value of their 
total number Nm = 27i'"''/Z„ to the total number of r- 
sites, n: 



7-)(™) 

am . 

n^oo nZn 



(3.7) 



similar to equation (2.5). These concentrations must sat- 
isfy the conservation law 



l = ^(fc-f l)a;fe. 



(3.8) 



k>0 



which states that the total number of compound spins, 
zeros and ones, is equal to the total number of r- 
sublattice sites. 

The set of equations (^)-( |3!^ ) allows us to perform 
a straightforward numerical analysis. However, let us 
try an analytical approach by assuming that |/int(fc) — 
/int(l)| ^ T for all fc > 1. Practically, this means that 
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T» |/i„t(l)|, 



(3.9) 



which is well satisfied at r>0.15, as we shall see in sec- 
tion IV. In this approximation one can set Wk = w for 
all k > 1. By subtracting Zn from 2„+i, we rewrite 
Eqs. ( |3.3| ) in a simple form (n > 1) 

Z„+i - (1 + 2wo)Z„ + {2wo - wi)Z„_i = 0. (3.10) 

Analogous equations for I?^"-* and are similar, but 
have different r.h.s. and boundary conditions: 



V, 



(0) 

n+l 



= 2wn{Zn - Z„_i), n > 1 



(3.11) 



(0) 



Vi% - (1 + 2^o)I?i'^ + (2«;o - wi)Vi:i, 
= wi(Z„_i - Z„_2), n>2 



0, 1?^°^ = 2wo. 



(3.12) 



V, 



(1) 



I? 



(1) 



0, 



As in the case of Eqs. (2^)- (2^), we look for solutions 



Zn, I?!"'' and Vn^ of the form 



c+\l + C-Xl, {a+ + b^°\i)Xl + (a_ + b'^>n)\'l, 
and (a+ + b^^^n)Xl + (a_ + b^}^n)X":, 

respectively. By inserting this ansatz we obtain 
X+=wa + 1/2 + y^2w2 + 1/4, 

and the concentrations of chain fragments of length k = 1 
and 0: 



Xi 



^(1) 



u;i(A+-l) 



A+(A^ + wi - 2it;o) ' 



Xq 



b'-^^ _ 2^0 A_| 



(3.13) 



(3.14) 



Generalization of Eg. ( 3.13 ) for to > 1 is straightfor- 
ward: The analogue of ( 3.12 ) now reads 



- (1 + 2t.o)I?i"' + {2wo - w^)vt\ 
= wi(Z„_„ - Z„_,„_i), n>TO + l, (3.15) 

while the boundary conditions can be written as 



(m)_ 



Wi. 



The abovementioned asymptotic behaviour of Z„ and 
I?!'"'' allows us to determine the concentration Xm as 



,(m) 



1 

(a7)^ 



■ Xl. 



(3.16) 



It is not difficult to check t he va lidity of th e sum rule 
(U) with equations ( |3.14| ), ( |3.13D and ( [3.161) . 

Just comparing two sets of equations, (3.12) and 



(1) 



( 3.15| ), one can conclude that the knowledge of all 
allows us to calculate any expression of higher rank, e.g. 



V, 



(m) 
n+m- 



-75(1) 



(3.17) 



At lower temperatures it is necessary to bring in more 
Boltzmann factors. The simplest extension of the tem- 
perature range 



T» |/int(2)| 



(3.18) 



does not impose restrictions on |/int(l)| any more, the 
latter is not necessarily much smaller than the tempera- 
ture. This leads us to the introduction of three Boltz- 
mann factors, wq, wi, a nd W k = u'2 ~ Wi exp[/int(l)/T'], 
if > 2 (see definition (3^)). Numerical results of the 
next section will show that such a description is a good 
approximation for T > 0.04. This is a systematic way 
to extend our approach to lower temperatures. By us- 
ing four Boltzmann factors, the system can be described 
down to T w 0.01. 

In Appendix A we outline how the method developed 
above works in the general situation, where temperature 
is restricted from below: T 3> |/int(fc)|- In this case we 
deal with a set of statistical weights 

{wo,wi,W2,. ■ . , Wfe_i, w„ ^Wk{n> k)}. (3.19) 



IV. NUMERICAL COMPUTATIONS 

A. Chain fragments: Zero-field results 

By using complete matrix diagonalization we have 
computed all energy eigenvalues of finite chain fragments 
up to fc = 6. From these eigenvalues, the exact val- 
ues of the free energy for various temperatures have 
been calculated. We decompose the free energy into 
an affine contribution fc/oo + /o and an "interaction" 
contribution fint{k). This decomposition is defined by 



lim /i„t(fc) =0. 



k 


T = 0.02 


r = 0.3 


T = 


1.0 





0.454050 


0.245974 




0.027743 




1 


-0.074381 


0.002104 




0.000148 




2 


-0.012040 


-9.540148- 


10-^ 


8.060468 


•10-^ 


3 


-0.000113 


-5.336452 • 


10-^ 


4.404270 


•10-9 


4 


0.000476 


-8.747215- 


10-8 


2.376055 


•10-" 


5 


0.0 


0.0 




0.0 




6 


0.0 


0.0 




0.0 





TABLE I. The "interaction" contribution /int(fc) of the free 
energy for T = 0.02, T = 0.3, and T = 1.0. 
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Table | shows computed values of /int(fc) for T = 0.02, 
T = 0.3, and T = 1.0. Free energy values for k — 5 and 
fc = 6 have been used to determine /o and /oo, so 

/int(5) = /i„t(6) = 

for all temperatures. The accuracy of this approximation 
can be estimated from figures ||, ||, and ^, where /int(fc) 
decays quickly and smoothly to zero. 



0.1 



0.05 



-0.05 



-0.1 




T=0.02 
T=0.1 



1 2 3 4 5 6 

k 

FIG. 2. Computed free energy vs k after subtracting the 
afRne contribution kfoo + fo for low temperatures T = 0.02 
(full) and T = 0.1 (dashed) 




FIG. 3. "Interaction part" of free energy as a function of k 
for high temperatures T = 1.0 (full) and T — 0.5 (dashed) 

For low temperatures (T < 0.3), the fc-dependence 
of the free energy is qualitatively similar to the ground 
state energy. This is illustrated in figure ||. As in 
the ground state energy, there is a clear minimum at 
k = 1, which corresponds to the phase transition se- 
quence (0) (1) (oo) at T = 0. 

For T > 0.3, temperature has smeared out the ground 
state structure completely. Deviations from the affine 
contribution kfoo + fo decay exponentially fast as a func- 



tion of fc. This is shown in figure I for T 0.5 and 
T = 1.0. 

In the crossover region T ~ 0.3 we observe an interest- 
ing phenomenon: As shown in figure ^, the minimum of 
the free energy has moved from fc = 1 to fc = 2. This is 
due to the spin-degeneracy of the chain fragment. The 
total spin of a chain fragment of length fc is Sp = (fc — 1)/2 
for low temperatures, so the degeneracy is 2sp + 1 = fc. 
This yields 



/do 



-T Ink 



as an additional contribution to the free energy. If we 
subtract this additional contribution from /int, the min- 
imum at fc = 1 is restored. This contribution is only 
relevant at intermediate temperatures: At low tempera- 
tures it is suppressed by the prefactor T, at high tem- 
peratures the In fc contribution to the spin entropy is no 
more dominant. Indeed, the magnetic behaviour of short 
chain fragments cannot be described by a single param- 
agnetic spin Sp. In other words, at high temperatures the 
spin-spin correlation length becomes smaller than fc. 



fin, 0.001 




0.0005 



-0.0005 



FIG. 4. "Interaction part" of free energy as a function of k 
at the crossover temperature Tc — 0.3 



B. Chain fragments: Magnetic susceptibility 

By including small magnetic fields into the Hamilto- 
nian, we have computed the zero-field susceptibility for 
fc = 1, . . . , 5 and various temperatures. 

There are two contributions to the magnetic suscepti- 
bility: 



1. The paramagnetic contribution, which is 



1 



(4.1) 



for low temperatures. 



2. The contribution Xa due to the antiferromagnetic 
correlations inside the chain fragment. Though the 
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exact form of this contribution is unknown, it is 
expected to be approximately linear in k and only 
slowly varying as a function of temperature. 

Figure ^ shows the magnetic susceptibility as a function 
of inverse temperature for fragment lengths k = 2, . . . ,5. 
The numerically computed values (various symbols) are 
almost perfectly connected by the exactly known para- 
magnetic contribution ( [4.1| ) (lines), where (fc — l)/2 has 
been inserted for Sp. Please note that there are no fit 
parameters. 




FIG. 5. Computed magnetic susceptibility f or fc — 2, ... ,5 
(symbols) and the paramagnetic contribution 
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FIG. 6. Computed magnetic susceptibil- 

ity for T — 0.02, . . . , 0.06 (symbols) and the paramagnetic 
contribution ( [4.l| ) (lines) 



atures can also be observed when plotting x as a function 
of the fragment length fc. This is shown in figure ^ for 
temperatures 0.02 < T < 0.06. 

For high temperatures, the correlation length becomes 
smaller than the fragment length. In this case, the effec- 
tive total spin can be interpreted as being composed of 
several independent blocks of individual spins. The typ- 
ical length of these blocks is the correlation length ^, so 
the effective total spin square Sp{sp + 1) is proportional to 
^fc, in contrast to the low temperature situation, where 
Sp = (k — l)/2. This results in a linear x vs fc depen- 
dence, which is numerically confirmed, as shown in figure 
^ There we have successfully fitted affine functions 



Xh = ao{T) + ai{T) fc 



to the computed data points. 

2.6 
2.4 
2.2 

2 
1.8 
1.6 
^ 1.4 
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0.8 
0.6 
0.4 

1 



(4.2) 







T=2.0 ■> 




T=1.5 + 




T=1.0 = 









2 3 4 5 

k 

7. Computed magnetic susceptibility for 

0,1.5,2.0 (symbols) and fitted affine functions (4.2) 



FIG. 

T = 1 
(lines) 



If we would extend figure ^ to large values of fc, we 
would certainly observe a crossover from parabolic to lin- 
ear fc-dependence at sufficiently large fragment length fc. 
The crossover region depends on temperature. 

We finish this section by considering intermediate tem- 
peratures, where the crossover behaviour can be observed 
at finite fragment length fc < 6. Figure || shows this 
intermediate temperature region. The crossover from 
quadratic to linear behaviour is governed by a logarith- 
mic correction 



Therefore, at least for low temperatur es, the suscepti- 
bility is well described by the Curie law ( 4T ) . For higher 
temperatures Xv ^^^^ strictly linear in . Addition- 
ally, with increasing temperature the antiferromagnetic 
contribution Xa decays slower than Xp, so Xa may be- 
come relevant for larger T. 

Of course, the coincidence of x with Xp for low temper- 



X. = 6o(T.)-t-6i(r,)fc + 62(^0 Infc, 



(4.3) 



which fits x(^) successfully for temperatures around T = 
0.3. This behaviour presumably is reminiscent of the 
crossover in the free energy, which takes place at the same 
temperature range. 
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Values for foo and wq , ■ • • , have been obtained from 
the numerical results outlined in subsection A. 



FIG. 8. Computed magnetic susceptibility for intermediate 
temperatures (symbols) and fitted functions (4.5) (lines) 



C. Mixed chains as ensembles of chain fragments 

This part of the paper is some kind of a synthesis of 
the analytic approach developed in Section III and nu- 
merical computations performed in Section IV. We shall 
illustrate this by computing the specific heat for two tem- 
peratures, low and intermediate, T = 0.02 and T = 0.3, 
and for reasons, which will be explained below, the mag- 
netic susceptibility versus j7o will be calculated for inter- 
mediate temperatures T = 0.1 and T = 0.3. Please note 
that in Section II we derived an analytical expression for 
the magnetic susceptibility, which holds in the vicinity of 
the (0) — > (1) transition, where all contributions of chain 
fragments of length k > 2 can be neglected, i.e. at low 
enough temperatures. 

The specific heat per compound spin is calculated as 



C 



-T 



where / is the total free energy of the alternating chain. 
In the thermodynamic limit, the partition function is 
given by the loga rithm of the maximum root Amax of the 
polynomial (A.3). The "reference state" in the derivation 
of the partition function is the perfect (oo)-structure, so 
the complete free energy is given by 

/ = -r ln(A,„ax)+/oo. (4.4) 
According to Table I we have the estimate 



|/int(fc)| < 0.02 for fc > 3. 

Therefore it is sufficient to work with four different Boltz- 
mann weights, wq, . . . ,W3. (At T = 0.3, this approxi- 
mation is even more accurate.) The polynomial (A.3) 
specializes to 

A"- (l + 2wo)A^+ {2wo~ wi)X^+ {wi - w^)\ 




FIG. 9. Specific heat of the full alternating chain as a func- 
tion of Jq sX T = 0.02. The insert shows the fine double 
peak structure around ^Jq^'""'. Note that for > —0.7 the 
specific heat approaches a non-zero value (temperature de- 
pendent!), which is the specific heat of the perfect alternating 
chain ... - 1/2 - 1 - 1/2 - 1 - .. .. 
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FIG. 10. Specific heat of the full alternating chain as a 
function of Jo at T = 0.3. It decays to zero for large negative 
Jo, e.g. C « 0.02 at Jo = -4. 

Figure || shows the resulting specific heat as a function 
of Jo for T = 0.02. As expected, the most interesting 
features can be found in the regions around Jq^'^"^ and 
j.{i,oo) ^ In order to understand the double-peak struc- 
ture, which is shown in the insert, consider the expression 




-(w2— Ws) = 0. 



/(T) = ^a;fc(T) ^^(r) 



k>0 
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for the total free energy of the alternating chain. As in- 
dicated, both the concentrations Xk and the free energies 
4>k of the finite fragments depend on temperature, so the 
specific heat involves Xk itself and its first and second 
derivatives. At low temperatures, Xk is a step-like func- 
tion, which qualitatively explains the splittings of the 
peaks and their shapes. 

At T = 0.3, temperature has smeared out the fine 
structure around Jq^'^^ and J'q^'°°'' . Both regions are 
not isolated from each other, nevertheless the double- 
peak structure is still visible in figure ^ 

Figures ^ and ^ illustrate the dependence of the spe- 
cific heat on the interior coupling constant j7o. The 
temperature-dependence of C will be published in future 
work. 

Calculation of the magnetic susceptibility per com- 
pound spin is straightforward. It is simply given by 



X ^ 



Xk, 

k>0 



(4.5) 



where Xfc is the susceptibility of a finite chain fragment of 
l ength k. The concentrations Xk according to ( [A.8| ) and 
( A. 10 ) are computed by specializing to only four differ- 
ent Boltzmann weights, xo is the susceptibility of an iso- 
lated paramagnetic spin-I/2. In units of /i^, xo = l/4r. 
Xi, . . . , X5 known numerically from subsection A. For 
Jo <_Jn^'°°'' the largest eigenvalue Amax ^ 1, so because 



of (A. 10) we can safely cut-off the summation (4.5) af- 
ter fc = 5. As Jo approaches J'q^'°°\ all concentrations 
Xk for finite k quickly decay to zero, as (oo) becomes 
the dominating configuration. The susceptibility Xoo of 
this configuration is unknown for T — 0.02, but we can 
extrapolate it easily for those temperatures at which x 
versus k is in the linear regime. From the calculations 
performed in subsection B, Xoo per compound spin can 
be estimated as 



T 
T ■ 



0.1 
0.3 



X' 

x' 



liuik- 
liuik- 



>oo ^ Xk 
>oo "kXk 



6.70, 
1.60. 



(4.6) 



Unfortunately, the convergence of the series ( |4.5| ) is very 

bad if Jo > Jo^'°°^ ■ At intermediate temperatures, this 
problem can be circumvented by making use of the sum 
rule (3^) and the relations 



Xkoc k 



1 for increasing Jo, 
for large k 
for large k. 



Amax 

Xk (X A,7^^x fo'' large k and 



This procedure is illustrated in Appendix B. However, in 
order to obtain reasonable results one needs to know the 
whole set of Xfe before the /c-dependence becomes linear. 
There is no problem in case of temperatures T = 0.3 and 
r = 0.1 (cf. figure 11), but the set of Xfc at T = 0.02 is 
too far from the linear regime (see figure |^) . 

In figures ^ and which show the magnetic suscep- 
tibility for T = 0.1 and T = 0.3, we have used figure || as 



a prerequisite, i.e. we extracted the fragment susceptibil- 
ities^, . . . , X5, and the asymptotic slope x' ■ Figures p] 
have the correct asymptotics: At large negative 
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and 

Jo the susceptibility approaches at large positive Jo 
it is in accordance with (4.6). 
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FIG. 11. Magnetic susceptibility of the full alternating 
chain as a function of ^Jo at T = 0.1. 




FIG. 12. Magnetic susceptibility of the full alternating 
chain as a function of at T = 0.3. 



V. SUMMARY AND OUTLOOK 

In this work we have combined analytical and numer- 
ical methods for calculating thermodynamical proper- 
ties of alternating Heisenberg chains. Interest in these 
spin systems is not purely academic. In fact, recent 
progress in the observation of quantum effects in quasi 
one-dimensional antiferromagnets, see e.g. ||^ and is 
commonly contributed by physicists and crystal growth 
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experts. It is out of our goal to recommend some definite 
compounds, which might be good candidates for mixed 
chains with such a competitive behaviour. Schematically, 
their simplest representation is given in our description 
of models A and B (see also illustrations in 

It is necessary to emphasize, that the unusual thermo- 
dynamic behaviour described in this paper can be effi- 
ciently observed in those compounds, where the values of 
the coupling constants Jo and Ji lead to competitive in- 
teractions. Then a major contribution to thermodynamic 
quantities comes from thermal changes of the fragment 
concentrations xj.. 

1. To make the basic ideas transparent, we employ a 
simple representation of these chains: "compound 
spins" and spin-1/2 sites alternate within the chain. 
They interact via a coupling constant J7i. A com- 
pound spin consists of two spin-1/2 sites forming a 
dumbbell configuration perpendicular to the chain 
direction. The constituents of a dumbbell interact 
via the internal coupling constant JTq. The intrinsic 
state of a compound spin is either spin-1 or spin-0, 
which gives rise to the hidden Ising symmetry of 
the Heisenberg chain. Spin-0 states can be treated 
as "defects" with respect to the perfect alternating 
chain ... — 1/2—1 — 1/2—1 — .... These equilibrium 
"defects" simply break the perfect chain into a set 
of non-interacting chain fragments of finite lengths, 
whose general formula is (i, 1)*^^ for length A: > 0. 

2. At zero temperature, one of three ground states is 
realized, depending on the value of Jq {Ji, the an- 
tiferromagnetic exchange, is set to —1): 

(0) , the regular structure of chain fragments of 
length fc = for j7o < — 2, this means that all com- 
pound spins are zero; 

(1) , the regular structure of chain fragments of 
length k — 1 for — 2 < Jo < —0.910, i.e. compound 
spins and 1 alternate with each other; 

(oo), only one "fragment" of infinite length for 
— 0.910 < Jo, i-e. all compound spins are 1. 
Although the first order transitions governed by Jo 
are smeared out at non-zero temperature, thermo- 
dynamic quantities display a strong dependence on 
Jo and temperature in the vicinities of the zero- 
temperature critical values. 

3. As a function of Jo, the specific heat displays 
two peaks, which are reminiscent of the zero- 
temperature transitions (0) (1) — > (oo). As 
shown in figure these peaks are split at low tem- 
peratures. The origin of this fine structure is the 
step-like behaviour of the fragment concentrations 
Xk , predominantly xq and xi . With increasing tem- 
perature, the peaks become broader and overlap, 
but the local minima still point out the vicinities of 
the critical Jo values, cf. figure |l^. The crucial de- 
pendence of the magnetic susceptibility on Jo and 
T is illustrated in figures |ll| and 12, Technically it 



is more difficult to compute yY Jo) because of the 
poor convergence of series (4.5). We have elimi- 
nated these difficulties by making use of the sum 
rule (U). 



4. All these pretransitional phenomena, which can 
be identified clearly, are provided by the config- 
urational contributions to thermodynamic quanti- 
ties. We have developed an analytic approach for 
calculating the free energy and the concentrations 
of chain fragments of different lengths. This is 
achieved by solving polynomial equations, whose 
degree increases at lower temperature. In fact, the 
hierarchical structure of the equations is regulated 
by /int(^), i-e. the part of the free energy of a chain 
fragment of length k which can be interpreted as 
the interaction of fragment ends. From numerical 
computations we observed that /int decays rapidly 
with k. 

5. For other alternating chains, where the non- 
compound sites are occupied by higher spins s 
instead of 1/2, there are several succesive zero- 
temperature transitions with increasing Jo- The 
last one is (4) (oo) if s > 3/2. For these chains, 
the occurence of several peaks in the specific heat 
would be possible. However, very low temperatures 
are needed to isolate some peaks: Although the suc- 
cessive transitions (0) — *■ (1) and (1) (2) are suf- 
ficiently separated, the distance between (3) (4) 
and (4) -> (oo) is only A Jo ^ 10"^. 

6. To investigate the thermodynamics of alternating 
chains with s > 1/2, other methods could be used 
instead of exact numerical diagonalization. A mod- 
ification of the spin- wave approximation, which can 
be adjusted to our problem, has been proposed by 
M. Takahashi Q. In Q and the quantum trans- 
fer matrix of exactly solvable spin chains has been 
investigated in order to calculate thermodynamic 
quantities. A combination of this approach and nu- 
merical methods could be applied to our problem. 

7. Generalization to higher dimensions is straight- 
forward: For instance in two dimensions, the 
"white" squares of a checkerboard lattice contain 
single spins, while the "black" ones are occupied 
by compound spins. Competition may again re- 
veal a hidden Ising-like variable and possible two- 
dimensional superstructures. Of course, similar 
models can be constructed on any bipartite lattice. 
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APPENDIX A 

As shown in section IV, /int(fc) decays rapidly at 
high temperatures. But even at very low temperatures, 
r = 0.02, we have the estimate 

|/int(fc)| « T for k> 3. 

According to the method developed in section III, it is 
sufficient to introduce the set of Boltzmann weights: 

{Wi, W2 = wi exp[(/iiit(l) - /int(2))/T], 
W3 = wi exp[/i„t(l)/r], Wn = wz{n > 3)}. 

In this appendix we outline the procedure for the gen- 
eral case, where all statistical weights (3.19) are taken 
into account. 

If we substitute n by n — 1 in ( |3.3| ) and subt ract the 
resulting expression from the original equation (3.3) we 
arrive at 

Zn+l- (1+ 2wo)Zn+ {2wq- Wi)Z„_i + (wi" W2)2„^2 
+ (W2- W3)Z„-3 H h {wk-1- Wk)Zn^k = 0, 



which is a generalization of (3.1C). It can be rewritten in 
compact form 



7?.fc+iZ„+i =0, n > fc, 
if the linear operator 



(A.l) 



"^fc+l/n+l = (l + 2lt;o)/n+ (2^0- W\)jn-\ 

+ {wi ~ W2)fn-2 H h (Wfc-1 - Wk)fn-k (A.2) 

is introduced. Directly related to the operator TZk+i is 
the characteristic equation 

Pk+i{\) = A'=+i - (1 + 2u;o)A*^ + (2u;o - u;i)A^-i 

+ (wi - u;2)A''"^ + • • • + {wk-i - Wfc) = (A.3) 

whose maximal root, A„iax, from the set {Ai, . . . , Afc+i} 
determines the asymptotic behaviour of the partition 
function. 

Analogous equations for the functions can be de- 
rived and written in compact form as 



Zri—i 



(A.4) 



where Wm must be substituted by 2wq for m = 0. 

In our formalism it is sufficient to find out the P^™) 
functions up to rank k (m = 0, 1, . . . , fc). For higher rank 
m>k one obtains 



(m) 



(A.5) 



This rela tionship becomes clear from the form of 
Eqs. (A.4) for m> k, where Wm — Wk- 



An obvious ansatz to solve Eqs. (A.l) is 



fc+i 

-^n — ^ ^ CiAj , 
1=1 



(A.6) 



whereas for solutions of Eqs. ( |A.4D we use 

k+l 

2?i''^=E(«r^+^r^")^"' r = 0,l,...,fc. (A.7) 

i=l 

We do not need the complete set of coefficients 
{a^l\b'f\ Ci} , which can be calculated by using the 
boundary conditions, but rather the set of ratios 
{6max/cmax}, which yiclds the concentrations 



pare the leading terms, whic 



After inserting the ansatz (A.7) into Eq. (|A.4D, we com- 



1 are given in terms of pow- 



ers of Amax- Neglecting all other terms is equivalent to 
taking the thermodynamic limit. The basic result for the 
concentrations reads 



r,(™) 



Wm(A„ 



1)A: 



fc— m— 1 
max 



where 



Qfe+i(A) = A 
-2{wi- W2)\''-^ 



max ) 

- k{wk-i - Wk) 



(A., 



^ dX 



Afc 



Combining Eq. 



with result (A. 8) yields 
1 



Xk+r. 



A" 



-Xk- 



(A.9) 



(A.IO) 



APPENDIX B 

For the temperatures of interest we can write the sus- 
ceptibility in the form 

X = xo + a;i xi + ^2 X2 + X3 (xa + X4 A^^x 

+X5 (Amix + KiIk + ■ ■ •) 

+X' A-L(l + 2A„L + 3A-4 + •••))• (B.l) 

The most singular behaviour is contained in the term 
which is proportional to x' ■ It can be accurately evalu- 
ated by u sing the sum rule (3.8). Let us reorganize the 
terms of (B.l) into three groups. The first one contains 
the non-singular contributions 



xo XO + X1X1+ X2 X2 + X3 (X3 + X4 A^^,, 
+x' A„iax(l + 2A„jax + 3Ajjj^^)). 

The second one is pseudo-singular: 

X3 X5 



1) 



(B.2) 



(B.3) 
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Although the denominator in (B.3) vanishes if A,nax ^ Ij 



the complete term remains finite, since 0:3 —> in this 
limit, too. (cf. Eq. (A. 8).) Due to the sum rule 



2:0 + 2xi + 30:2 + (4 + 5Aj^j^^ + 6A^ax + • • 
the last contribution 

can be rewritten in the form 

x' Amax(l -xq- 2xi - 3x2), 
which evidently approaches x' when Amax ~* 1- 



■) = 1, 



(B.4) 
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